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Introduction 

M ANY computational fluid dynamics codes for turbo- 
machinery use the Baldwin -Lomax (B-L) turbulence 
model. 1 It is easy to implement in two dimensions and works 
well for predicting overall turbomachinery performance. How- 
ever, it is awkward to implement in three dimensions, often 
has difficulty finding the length scale, has a crude transition 
model, and neglects freestream turbulence, surface roughness, 
and mass injection. 

The k-co model developed by Wilcox 2 is an appealing alter- 
native for several reasons. First, it is the only two-equation 
model that can be integrated to the wall without requiring 
damping functions or the distance to the wall, and hence, 
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should behave well numerically. Second, the effects of free- 
stream turbulence, surface roughness, and mass injection are 
easily included in the model. Finally, transition can be simu- 
lated using the low Reynolds number version of the model. 3 

Menter applied the k-co model to external flows and showed 
very good results for flows with adverse pressure gradients. 4 
Liu and Zheng 5 described their implementation of the k-co 
model in a cascade code that included an area change term to 
account for endwall convergence. They validated the model 
for a flat plate, and compared the B-L and k-co models to 
measured surface pressures for a low-pressure turbine cascade. 
Since they did not use the low Reynolds number version of 
the model, their results showed problems resulting from early 
transition. 

In this Note the low Reynolds number k-co model was in- 
corporated in the author’s quasi-three-dimensional turboma- 
chinery analysis code. 6 The code includes the effects of rota- 
tion, radius change, and stream-surface thickness variation, and 
also includes the B-L turbulence model. The k-co model was 
implemented using many of Menter’ s 4 recommendations and 
an implicit approximate -factorization scheme described by 
Baldwin and Barth. 7 The model was tested for a transonic com- 
pressor with rotation and variable stream-surface radius and 
height, and for a transonic turbine vane with transition and 
heat transfer. Results were compared to the B-L model and 
to experimental data. 

Quasi-Three-Dimensional Navier- Stokes Code 

The quasi-three-dimensional Navier- Stokes equations have 
been developed for an axisymmetric stream surface in an (m, 
6) coordinate system as shown in Fig. 1. The meridional co- 
ordinate m is given by dm 2 = d z 2 + dr 2 , and the tangential 
coordinate 0 rotates with the blade row with angular velocity 
0. The radius r and the thickness h of the stream surface are 
assumed to be known functions of m . The equations were 
mapped to a body-fitted coordinate system, nondimensional- 
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Fig. 1 Quasi-three-dimensional stream surface for a compressor 
rotor. 

ized, and simplified using the thin-layer approximation. The 
flow equations and details of the explicit Runge-Kutta scheme 
are given in Ref. 6. 

k-eo Turbulence Model 

The k-co turbulence model solves two transport equations for 
the turbulent kinetic energy k and the specific dissipation rate 
<o. A basic formulation is for fully turbulent flows and a low 
Reynolds number formulation is used for transition. 3 The 
quasi-three-dimensional form of the low Reynolds number 
k-co equations is as follows: 


d,q + Ud ( q + Vd v q - G = -(/>- D) (1) 
pRe p 

q = [K o>) T (2) 

fi T = a*(pk/co) (3) 

The r and h appear in the 6 metrics and the Jacobian, and 
0 appears in the relative contravariant velocities U and V : 
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The diffusive terms G use the thin-layer approximation: 
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The production terms P use the vorticity Cl as suggested by 
Menter 4 : 
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The destruction terms D are given by 
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The baseline k-co model has five coefficients, /3 = 3/40, 


P* = 9/100, o- = 1/2, cr* = 1/2, and a = 5/9, and the trivial 
constant a* = 1. The low Reynolds number model replaces 
three of the constants with ft* = (9/100 )F <3 , a = (5/9 )(F a /F^), 
and a* = F^, where F, are the following bilinear functions of 
the turbulence Reynolds number Re T : 
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with q 0 = 1/10, a$ = p/3 = 1/40, and R p = 8 , = 27/10, 

F* = 6. 

Boundary Conditions 

At the inlet the turbulence intensity Tu and turbulent vis- 
cosity pi T were specified. Then co was found from Eq. (3) and 
k was found from 


k = iTu 2 Ul 


( 10 ) 


On solid walls k = 0, and co was set using Wilcox’s rough- 
ness model: 
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The equivalent sand grain roughness height k r was set to five, 
giving a hydraulically smooth surface. 

Values of k and co were extrapolated at the exit and treated 
as periodic across trailing-edge wake cut lines and between 
blade rows. 


Approximate Factorization Solution Scheme 

An implicit approximate -factorization scheme, described by 
Baldwin and Baith, 7 was used to solve the k-co equations. Ad- 
vective terms were approximated using first-order upwind dif- 
ferences and diffusive terms were approximated by second- 
order central differences. The destruction terms were linearized 
as suggested by Menter 4 and were included in the streamwise 
£ factor. The production terms were treated explicitly. The k 
and co equations were solved uncoupled from each other and 
from the flow equations. 

The CPU time for one iteration of the k-co model was only 
5% more than that for one application of the B-L model. 
Experience has shown that it is sufficient to update the B-L 
model every five iterations. In general, it was necessary to 
update the k-co model every two iterations with twice the time 
step of the flow solution, making a k-co solution about 18% 
slower than a B-L solution. 


Transonic Compressor Rotor Wake 

The transonic compressor rotor studied by Suder et al. 8 was 
used to test the quasi-three-dimensional effects in the model. 
This rotor was used for a three-dimensional blind test case for 
turbomachinery codes organized by ASME and the Interna- 
tional Gas Turbine Institute. Many of the codes used for the 
test case were unable to predict the measured wake spreading 
and decay. It was conjectured that the B-L model was re- 



1178 


J. PROPULSION, VOL. 12, NO. 6: TECHNICAL NOTES 



Fig. 2 Computed and measured near-wake Mach number pro- 
files for a transonic compressor rotor. 



Fig. 3 Computed and measured distributions of isentropic Mach 
number for the von Karman Institute turbine vane. 


sponsible for the poor predictions and hoped that the k-to 
model would improve the results. 

A section of the rotor at 70% span was analyzed. The 
stream-surface r and h were specified to model the converging 
endwalls. The radius was specified as a line at 70% of local 
span, and h was specified as the local span normalized by the 
inlet span. A meridional view of the stream surface is shown 
in Fig. 1. The choice of r and h can have a first-order effect 
on the solution. However, the algebraic guess used here often 
gives reasonable results that can be used to compare turbulence 
models to each other. 

A C-type grid with 319 X 45 points and a wall spacing of 
y + < 2 was used. The calculations were run 2000 iterations, 
which took about 3.25 min for the k-to model on a Cray C-90 
computer. 

Figure 2 compares computed and measured Mach number 
profiles 0.28 chords downstream of the trailing edge. The k-to 
(solid line) and B-L (dashed) results are nearly identical. The 
computed Mach number in the core flow is slightly high as a 
result of the specified stream surface, and the computed wakes 
are narrower and deeper than the measured wakes. Varying to ia 
by five orders of magnitude and varying Tu from 3 to 6% had 
little effect on the computed wake spreading. It is now thought 
that the measured wake spreading may be caused by unsteady 
vortex shedding from the trailing edge. 

Ttonsonic Turbine Cascade 

A transonic turbine vane tested by Arts et al. 9 was also com- 
puted. Comparisons were made with measured surface pres- 
sures and heat transfer data. A C-type grid with 383 X 49 
points and a wall spacing of y + < 1.5 were used. Surface heat 
transfer converged to plotting accuracy in 5000 iterations, 
which took about 8 min on the Cray C-90. 

Computed distributions of isentropic surface Mach number 
are compared to experimental data for exit Mach numbers of 
0.875 and 1.02 in Fig. 3. The B-L and k-to models gave iden- 
tical results. The subsonic results agree very well with the 
experimental data. The transonic results slightly underpredict 
the Mach number on the rear (uncovered) part of the suction 
surface. 

Figure 4 compares measured and computed values of sur- 
face heat transfer coefficient H [W/(m 2 K)] vs distance 5 (mm) 
along the vane surface for a case with exit M ex)l = 0.875, Re exil 
= 1 X 10 6 , and Tu = 4%. The experimental data (triangles) 
show laminar flow on the entire pressure surface, and on the 



Fig. 4 Surface heat transfer coefficient predicted by three tur- 
bulence models. 


suction surface up to transition at 5 = 60 mm. The laminar 
regions have augmented heat transfer caused by the high ffec- 
stream turbulence. This augmentation is not modeled by any 
of the turbulence models. 

The baseline k-co predicted fully turbulent flow on the suc- 
tion surface and showed transition near S = —20 on the pres- 
sure surface, giving high values of H. The B-L solution 
agreed closely with the data using the simple transition model 
proposed in Ref. 1. The low Reynolds number k-to solution 
also agreed closely with the data; however, the transition point 
was set by the choice of the inlet values of a>, as characterized 
by the inlet turbulent viscosity. The value of /x rin was varied 
by about three orders of magnitude to produce solutions that 
ranged from fully laminar to almost fully turbulent. Only a 
small range of values of ii TSa gave transition near the measured 
location. The results shown here used [/x r /jtL*J in = 7 X 1CT 4 . 
It is emphasized that although the location of transition de- 
pended strongly on fi TAtt or ti> jn , the heat transfer in the fully 
turbulent region was nearly independent of (u in . 

Concluding Remarks 

Wilcox’s k-to turbulence model has been added to a quasi- 
three-dimensional Navier- Stokes code for turbomachineiy. 
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The model included the effects of rotation, radius change, 
and stream-surface convergence. An upwind-implicit approx- 
imate-factorization scheme was used to solve the turbulence 
model equations uncoupled from the flow equations. The nu- 
merical scheme was robust, but about 18 % slower than the 
B-L model. 

Calculations were made of a transonic compressor rotor with 
significant quasi-three-dimensional effects. The results showed 
very close agreement between the B-L and k-to models, but 
both models failed to capture the measured wake spreading. 
Calculations were also made of a transonic turbine vane with 
transition and heat transfer. Both turbulence models showed 
very good agreement with measured surface pressures. Surface 
heat transfer was predicted reasonably well by the B-L model, 
considering the simple transition model used. The low Reyn- 
olds number k-co model gave similar results, but the predicted 
transition location was sensitive to inlet values of a>. Overall 
the k-<o model behaved well numerically, but predictions were 
not decisively better than those made with the B-L model. 
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